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Abstract 

In general, Picard-Fuchs systems in N = 2 supersymmetric Yang-Mills theories 
are realized as a set of simultaneous partial differential equations. However, if the 
QCD scale parameter is used as unique independent variable instead of moduli, 
the resulting Picard-Fuchs systems are represented by a single ordinary differen- 
tial equation (ODE) whose order coincides with the total number of independent 
periods. This paper discusses some properties of these Picard-Fuchs ODEs. In con- 
trast with the usual Picard-Fuchs systems written in terms of moduli derivatives, 
there exists a Wronskian for this ordinary differential system and this Wronskian 
produces a new relation among periods, moduli and QCD scale parameter, which 
in the case of SU(2) is reminiscent of scaling relation of prepotential. On the other 
hand, in the case of the SU(3) theory, there are two kinds of ordinary differential 
equations, one of which is the equation directly constructed from periods and the 
other is derived from the SU(3) Picard-Fuchs equations in moduli derivatives iden- 
tified with Appell's F± hypergeometric system, i.e., Burchnall's fifth order ordinary 
differential equation published in 1942. It is shown that four of the five indepen- 
dent solutions to the latter equation actually correspond to the four periods in the 
SU(3) gauge theory and the closed form of the remaining one is established by the 
SU(3) Picard-Fuchs ODE. The formula for this fifth solution is a new one. 
PACS: 11.15.Tk, 12.60.Jv, 02.30.Gp, 02.30.Hq, 02.30.Jr. 



I. INTRODUCTION 



It has been recognized that the low energy effective action of N = 2 supersymmetric Yang-Mills 
theory for any Lie gauge group including at most two derivatives and four fermions is dominated 
by a holomorphic function called prepotential Perturbatively, this prepotential is a sum of 
classical part and one-loop contribution, and further contributions from higher loop diagrams are 
excluded by the non-renormalization theorem in N = 2 theory. However, T was expected to be 
affected by instantons and hence it's non-perturbative determination was a longstanding problem. 

In the case of SU(2) gauge theory, i.e., without any quark hyper multiplet, Seiberg and Wittenlii 
showed that the vacuum configuration of the N = 2 action was parameterized by the moduli 
u = (ticf) 2 ) (<ft is a complex scalar field in the adjoint representation of the gauge group) and 
singularities on this parameter space (moduli space) could split into pieces by instanton effect. 
This instanton corrected moduli space is often called quantum moduli space and Seiberg and 
Wittenii identified the quantum moduli space with the moduli space of a certain elliptic curve of 
genus one. According to their ansatz, since the vacuum expectation value a = ((f)) and its magnetic 
dual can be also regarded as periods of a meromorphic 1-form on the elliptic Riemann surface, 
these periods can be calculated as a linear combination of solutions to Picard-Fuchs equations 
(for a historical review and introduction of Picard-Fuchs equation in mathematics, see Gray&). 
Once the periods are calculated, it is immediate to obtain the prepotential because of rigid special 
geometry. In this way, Klemm et aM determined the SU(2) prepotential. 

Seiberg and Witten's approach to N = 2 supersymmetric SU(2) Yang-Mills theoryiH 
was extended to other higher rank gauge group cases coupled with or without quark 
hypermultiplets,i'i0IS@iniJliJll0'E10 and it was found that the quantum moduli spaces of those 
gauge theories could be identified with those of Riemann surfaces with certain genus. In these 
studies, general algorithms to get Picard-Fuchs equations were develop ed,0Jll'lli'@il but these 
equations are in general realized as a set of simultaneous partial differential equations (PDEs) 
in terms of moduli derivatives. For this reason, it is not easy to solve Picard-Fuchs equations, 
especially, in higher rank gauge goup cases. 

However, if the QCD scale parameter instead of moduli is used as unique independent variable, 
the resulting Picard-Fuchs systems will be represented by a single ordinary differential equation 
(ODE) whose order coincides with the total number of independent periods. Then the problem 



solving Picard-Fuchs equations can be encoded into the language of ODE and therefore the study 
of periods are simplified. As another feature of this formalism, we remark that in contrast with the 
usual Picard-Fuchs systems written in terms of moduli derivatives there exists a Wronskian for these 
ordinary differential systems and the Wronskian produces a new relation among periods, moduli 
and QCD scale parameter. Especially, in the case of SU(2) it is quite reminiscent of scaling relation 
of the prepotential.0S0 This relation is highly non-linear in the case of higher rank gauge group, 
but it reflects the structure of the Picard-Fuchs ODE. Sec. II discusses these Picard-Fuchs ordinary 
differential systems. This realization of Picard-Fuchs systems via ODE becomes interesting when 
we consider a relation to hypergeometric differential equations in multiple variables. For example, 
in the case of SU(3) gauge theory, we can find another ODE which gives equivalent periods. That is 
Burchnall's fifth order equation directly constructed from Appell's hypergeometric differential 
equations.il In general, the dimension of the solution space of a single ODE constructed from 
simultaneous PDEs can exceed that of the original pdeH (see also Srivastava and KarlssonH and 
references therein), and Burchnall's equation is the case. In addition, the extra solution which is 
not a solution to the original PDEs is known to have a very characteristic form. In the case of 
Burchnall's equation for the SU(3) gauge theory, since four of the five independent solutions are 
found to correspond to the four periods of the SU(3) gauge theory (this identification is explicitly 
checked at the semi-classical regime) and these four periods are also solutions to the SU(3) Picard- 
Fuchs ODE, it is possible to extract a differential equation only for the fifth solution (although the 
fifth solution is irrelevant to the underlying physics). The formula for this fifth solution obtained 
in this way is a new one and takes a very different form compared with other fifth order ODEs 
constructed from a set of PDEs, e.g., a product of two Bessel functions or Whittaker functions (the 
"fifth solution" to these two cases are quite reminiscent each other). In our derivation, it is crucial 
to notice that the fifth solution of the fifth order ODE (that is, Burchnall's equation) associated 
with Appell's F 4 is constructed by "subtracting" fourth order equation satisfied by only periods of 
the SU(3) Seiberg-Witten curve which are part of the solutions to the fifth order ODE. In Sec. Ill, 
we discuss these aspects of Burchnall's equation as application of the SU(3) Picard-Fuchs ODE to 
a theory of hypergeometric equations. Sec. IV is a brief summary. 

Remark: When we simply say as "Picard-Fuchs ODE" , it always means a single ordinary 
differential equation in terms of QCD scale parameter derivatives. 



II. PICARD-FUCHS ORDINARY DIFFERENTIAL SYSTEMS 



A. The hyperelliptic curve 



Firstly, let us recall that the exact solution to the SU (n + 1) (n G N) gauge theory as an 
example. On the affine local coordinates x,y G C, the hyperelliptic curve and the Seiberg- Witten 
differential are given 



y 2 = w : 
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where z = A^J2 1} and 
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(2.2) 



Eq.( |2.2| ) shows the simple singularity with moduli Sj. This hyperelliptic curve can be compactified 
to a Riemann surface of genus n after addition of infinity, hence there must be non-contractible 
2n-cycles on this surface and these cycles can be chosen as the canonical bases, i.e., on fl <x,- = 
Pi n j3j = 0, cti fl /5j = —j3j n«j = +<5jj- Then we introduce the period vector 



n 



(2-3) 



where 



f A 5H /, a D . = f 



(2-4) 



Below, we often denote moduli as u = s 2 and t> = s 3 . 



B. Derivation of Picard-Fuchs ODE 

One of physically interesting behavior of periods is in the weak coupling region. Of course, 
the study of periods can be proceeded by using Picard-Fuchs equations, and in general periods 
are represented by series in moduli. However, by some rearrangement they are found to have a 
striking feature in their form. Namely, they can be summarized as 

periods = classical part + instanton corrections. (2-5) 



Eq. (|2.5| ) suggests that it is more convenient to construct periods as a series in QCD scale parameter 
A rather than moduli. In a sense, ( |2.5| ) might be implied by Ito and SasakuraS in their observation 
of general form of (a certain type of) Picard-Fuchs operators, which is summarized schematically 



as 



L = L d + AL A . (2.6) 

This equality means that the Picard-Fuchs operator L is a sum of operator L c \ whose kernel is 
classical periods and some operator La- Once the classical periods are known, we can calculate in- 
stanton corrected periods, that is, the kernel of L, by assuming La as a perturbation term for small 
A (at semi-classical regime). However, since in the case of other gauge theories with higher rank 
gauge groups Picard-Fuchs system is represented by a set of PDEs, such perturbative calculation 
involves technical problems, therefore another method to obtain instanton corrected periods should 
be developed. One of the candidates in view of differential equation is to construct Picard-Fuchs 
equations by regarding A as unique independent variable instead of moduli. Then the resulting 
Picard-Fuchs equations will be expressed by an ordinary differential equation. Since the QCD 
scale parameter always appears in any gauge theory with or without (massive) hypermultiplets, 
this formulation is convenient when we generalize the method to various gauge theories with any 
rank gauge group. 

Now, let us consider the derivation of this ordinary differential equation. In general, fc-times 
differentiation of Xsw over z gives 

d k X sw polynomial in x 

= ^Ti **> ( 2 - 7 ) 

but the right hand side can be decomposed into a sum of Abelian differentials and a total derivative 
term, if the well-known reduction algorithm is used.00S0ill Accordingly, collecting Q2.7D for 
various k can generate a differential equation. In addition, since the all independent periods should 
be solutions to this equation, the order of the equation must coincide with the total number of them 
and in fact it is determined as 2n. This reduction method is easily confirmed, if the Seiberg-Witten 
curves are hyperelliptic type. In this way, we get the Picard-Fuchs ODE in the form 



d 2n d 2n-l 



n = o, (2.i 



where functions in moduli. 



Also when massive quarks are included, we can obtain a similar ordinary differential equation, 
but in this case mass dependent polynomial appears in the denominator of the right hand side of 
( |2.7| ). Reduction of such massive differential was also recognized by Marshakov et a/.H in their 
construction of massive WDVV equations. 



C. Examples of Picard-Fuchs ODE 



Let us see examples of Picard-Fuchs ODE. The first one is the SU(2) case and then the coeffi- 
cients in (|2.8| ) are given by 
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Next, let us consider the SU(3) case. In this case, the coefficients are 
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where AsutS) is the product 



with the discriminant 
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(2.12) 



of the SU(3) hyperelliptic curve. 

We can easily obtain Picard-Fuchs ODE for SO(5) or Sp(4) gauge group in a similar manner. 
Though the Sp(4) hyperelliptic curve may be seen to be different from that of the SO (5) theory, 
the isomorphism of Picard-Fuchs equations between these two theories could be observed by Ito 
and SasakmaS Also in the case of Pieard-Fnehs ODE, this isomorphism can be easily established 
by the same transformation]^] 



D. Wronskian 



As is well-known, the scaling relation of SU(2) prepotential can be generated from Wronskian of 
Picard-Fuchs equation, but this is valid only for this SU(2) theory because Picard-Fuchs equations 
in other gauge theories consist of partial differential equations. In such theories, "Wronskian" does 



not generally exist. However, our Picard-Fuchs ODE (|2.8|) admits a Wronskian given by 
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where ' = d/dz. Substituting fl2.13|) into (|2.8|) shows that Wsu{n+i) satisfies 



W SU(n+l) + C 2n-\W S U{n+l) = 0, 



(2.14) 



which is integrated to give 
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where const, is the integration constant to be determined from the comparison of left and right 
hand sides of ( [2.1 5| ) by asymptotic behavior of periods, but may depend on moduli because in our 
formulation moduli are regarded as constant. 

Note that ( 2.15Q produces a new non-linear relation between periods and other parameters. 
For example, in the case of the SU(2) theory with the normalization used by Klemm et a/.,i (ETTH) 
gives 



/a 



(2.16) 



su{2) ~ ~ w 

which is a quite reminiscent expression with the homogeneity relation of prepotentiaLiSil Simi- 
larly, for the SU(3) theory, we have 
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where the integration constant is normalized to 1 for convenience in the next section. Note that 
the regular singularities of ( |2.17| ) are the same with those of SU(3) Picard-Fuchs ODE. 



III. BURCHNALL'S EQUATION 



A. Multi-term differential equation 

Picard-Fuchs equations obtained in the previous section can be shown to be classified in terms 
of multi-term ordinary differential equation discussed by Burchnall in his study of the relationship 
among hypergeometric differential equations in multiple variables and certain type of ordinary 
differential equations.il Here, multi-term ordinary differential equation is defined by: 

Definition : k-term ordinary differential equation is the differential equation taking the form 

i 



n = o (3.1) 



1=1 

for some I, where f and gi are polynomial differential operators in the Euler derivative 9 Z = zd/dz. 
k is the total number of f and non-zero gi. 

For example, in terms of the Euler derivative, our SU(3) Picard-Fuchs ODE can be rewritten 
as the four-term ODE with 



f = 2916(x-y)\-l + 9 z ) 2 9 z 2 , 



g l = -81 (x - y) 9/ 43x - 52y - 60 (x - y) 9 Z + 4 (23x + 13y) 9, 
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g 2 = g [(386x - 251?/) 9 Z 2 + 144 (2x - ly) 9 Z 3 + 36 (19a; + y) 9 Z 4 - (x - y) (10 + I39 z 
g 3 = -5 (1 + 39 z ) (2 + 39 z ) (-1 + 60,) (1 + Q9 Z ) , (3.2) 

where x = 4u 3 /27 and y = v 2 . 

According to Srivastava and Saran^ who extended the work of Burchnall to four-term ODE, 
our SU(3) Picard-Fuchs ODE seems to be representable by a hypergeometric function in the 
homogeneity form F(pz,qz,rz), where p,q and r are parameters. In this paper, we could not 
specify this function, but since the kernel of the SU(3) Picard-Fuchs ODE is essentially written 
by Appell's _F 4 function, some property of F 4 may appear in our SU(3) Picard-Fuchs ODE as 
a four-term equation. Furthermore, more detailed study indicates that Picard-Fuchs ODE in 
any rank gauge group can be classified as fc-term equation, but k seems to correspond to 2 x 
(rank of the gauge group). 

Finally, note that the Picard-Fuchs ODE in SU(3) gauge theory has a factor (—1 + 9 Z ) 2 9 2 in 
/-polynomial. Therefore, the indicial indices at semi-classical regime are degenerated to —1 and 0. 



This indicates that there are logarithmic solutions at this regime. Of course, similar observation 
holds also for SO(5) and Sp(4) Picard-Fuchs ODEs. 



B. Appell's equations and Burchnall's equation 

As is well-known, a hypergeometric function admits a lot of transformations and reducibilities. 
For example, Gaussian 2 F\ system has 24 solutions and Appell's F\ has 60 solutions. However, 
these solutions can be more systematically constructed, if we consider an equivalent ODE. In fact, 
Srivastava and Saran succeeded to find 120 solutions to the ODE for Lauricella's function.il 
Also in this sense, study of ODE for hypergeometric partial differential system is interesting. The 
method used in these studies followed to Burchnall's work.il In this paper, we do not attempt to 
obtain all solutions solutions to Burchnall's equation (see below) like Kummer's 24 solutions, but 
we can show that the basic five solutions to Burchnall's equation, especially, the extra solution 
which is not a solution to system, can be derived by using SU(3) Picard-Fuchs ODE. Of course, 
it may be interesting if also this extra solution can be represented by Appell's F4, but we do not 
know whether it is possible or not. Nevertheless, we can establish the fifth solution as a simple 
formula. The reader should notice that our method presented in this paper is to use a fourth order 
ODE (SU(3) Picard-Fuchs ODE) satisfied by periods of Riemann surface in genus two (SU(3) 
Seiberg-Witten curve) and therefore our method is quite different from those mentioned above. 

Firstly, let us recall that the SU(3) Picard-Fuchs system! 



n = 0, 

n = 0, (3.3) 



27 V y 6 J V ' y 6, 
where we have introduced x = 4m 3 /(27A| [/ ^ 3 ^) and y = v 2 /Ag U ^, and 9 X = xd/dx and 9 y = yd/dy 
are Euler partial derivatives. In the terminology of previous subsection, ( |3.3j ) consists of two two- 
term equations and is nothing but the Appell's differential system for the type _F 4 hypergeometric 
double series 

fl(<.,ft7,7 ;*,») = JE (7)m(7% ^ (3.4) 

where a = (3 = — 1/6, 7 = 2/3 and 7' = 1/2. However, by the scaling transformation 



we see that (13.3) turns to 



[0 X (0 X + j-l)-xz(B x + 9 y + a){9 x + 9 y + p)]F = 0, 

9y(9 y + l' - 1) - yz(0 x + 9 y + a){9 x + 9 y + (3)] F = 0, (3.6) 
whose one of analytic solutions near (x, y) = (0, 0) is given by 

At first sight, since F(xz, yz) reduces to _F 4 for z — > 1, this scale transformation may be trivial, but 
( |3.6|) was used as a starting point in the Burchnall's work on a set of partial differential equations.il 
In fact, Burchnall noticed on the homogeneity relation of F 

(0 X + 6 y - 9~)F = 0, (3.8) 

where 9~ is the ordinary differential operator 9~ = zd/dz, and finally arrived at the ordinary 
differential equation of fifth order (see also appendix A) 

/o - 2(x + y)zf 1 (e 7 +a)(e~ + (3) + l -{x - y)zf 2 (9~+a)(9j + P) 
+(x - yfz 2 h{9~+a){9~+a + /?)(% + + 1)]F = 0, (3.9) 

where 

/ O = 0~(0~ + T -1) (e 7 + 1 '-l) (9~ + 1 + 1 '-2) (^+| + y-2), 

/2 = ( 7 _V)( 7 + y_2)( %+ 2 + r_i), 

In contrast with the SU(3) Picard-Fuchs ODE, since Burchnall's equation ( |3.9| ) is classified as a 
three-term ODE, it reflects a consequence of the nature of two variables hypergeometric function 
in the form F(pz, qz), where p and q are parameters.il 

As a direct calculation shows, ( ^.9| ) can not be expressed in the form [9~+ K(x,y, z))LF = 0, 
where K(x, y,z) is some function of x, y and z, and L is some differential operator of fourth order. 
Accordingly, such L does not exist and therefore the relation between our SU(3) Picard-Fuchs 



ODE and Burchnall's equation is non-trivial. However, note that the coefficient of highest power 
in the Euler derivative corresponds to the discriminant of the SU(3) curve 
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Therefore, in a sense Burchnall's equation is reminiscent of the SU(3) Picard-Fuchs ODE, but 
these two are not completely equivalent. Clarifying the relation between these two equtions is the 
subject in the rest of the paper. 



C. Four solutions at semi-classical regime 

It would be instructive to get solutions explicitly around A S u(3) — (z — oo) which is a regular 
singular point of the equation. In the work of Burchnall, solutions of ( |3.9|) were not calculated at 
any singularities, but since ( |3.9| ) is a linear ordinary differential equation, it is easy to solve it by 
traditional Frobenius's method under the assumption F = z~ v J2^=o A n z~ n for some v and A n . 
Then the indicial indices are determined as 

v = a, /3, a + 1, (3 + 1, + +2) (3.12) 

or equivalent ly, 

1 5 19 . 

vx = -g, »2 = g^3 = ^ (3.13) 

where v\ and z/ 2 are actually double roots. The solution for z/ 3 is the subject in the next subsection. 
Eq. flOD produces the recursion relations 

(x - y) 2 piA 1 - fa - a){vi - (3)(j)iA = 0, 

(x - y) 2 p n A n - {vi + n- oi-l){Pi + n- ft - l)cr n A n _i + XnA n - 2 = 0, n > 1, (3.14) 

where 

p n = (2ui + 2n - j — — 2)(vi + n — a)(vi + n — (3)(vi + n — a - l)(vi + n - (3 — 1), 
a n = (ui + n-'y - 7 )(2i/» + 2n - 7 - j')[x(ui + n - 7 - 1) + y{vi + 72-7-1)] 

+(ui + n- l)(2vi + 2n - 7 - 7' - 2)[x(vi + n - 7) + y{v { + 72-7')], 
Xn = {n + n - 2)(ui + n - 7 - + n - 7' - l)(vi + n - 7 - 7 ){2vi + 272-7-7') (3.15) 



with Aq = 1. Here, repeated indices are assumed not to be summed. If these recursion relations 
are used, the solutions corresponding to respective indicial indices will be obtained, but we must 
be careful, because there are indicial indices which differ by unit among them, i.e., a and a + 1, 
and (3 and /3 + 1. For example, let v = a. Then the recursion relations produce the n-th coefficient 
as a linear combination of Aq and A%. Thus the solution is given by a linear combination of 
z~ a (A + • ■ ■) and z~ a ~ 1 (Ai + •■•), but the "indicial index" of the last series can be seen as a + 1. 
Therefore, the last series can be also regarded as a solution corresponding to this index. In fact, 
it is easy to see that explicit construction of the solution supports this observation. Of course, in 
this case, since we would like to get a solution corresponding to the index a, A\ can be set to zero 
without loss of generality. For this reason, A\ is chosen as zero for the indices a and (3, while that 
for a +1,(3+1 and (7 + 7' — 2)/2 should be determined from the first equation in (|3.14 ). 
In this way, we get the regular series solutions (i = 1, 2) 



00 
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where the first few coefficients are given by 
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On the other hand, the degeneracy of v\ and z/ 2 produce the logarithmic solutions (j = 1, 2) 

(Pi = (fj In - + z~ v i B hn z~ n , (3.18) 
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where some of Bj n are 
Bia = 0, 
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It is interesting notice that these series are composed by a series in powers of l/{x — y) = 27/(4ti 3 — 
27v 2 ) which detects the discriminant of the semi-classical SU(3) curve. This feature is useful when 
we consider the structure of the quantum moduli space of the SU(3) gauge theory. 

For this purpose, let us recall the work of Klemm et a/J In the course of the analysis of the 
quantum moduli space of the SU(3) gauge theory, Klemm et a/.@ found that the quantum moduli 
space could be better understood as the complex projective space CP 2 with singularities which 
correspond to the strong coupling regime. Then this space can be covered by the three local 
(inhomogeneous) coordinates 

p .(^ _v 2 \ p fiu\ A«\ . / 27,\27A*\ 

For this reason, periods should be obtained at each coordinate patch, hence the periods derived 
in this way are locally valid. However, our solutions have a (slightly) nice property, because the 
basis of solution space are common both on P2 and P3. That is, to get periods on P2, it is enough 
to further expand (fi by v , while on P 3 by large u. 

In fact, we can see that the four periods are expressed by linear combinations of ( |3.16| ) and 
( 13.181 ). For example, to match periods on the patch P3, let us define uji and Qj by linear combina- 
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Here, t/>(x) = c?lnr(a;)/da; is the digamma function and 2 -Pi is the hypergeometric function whose 
series representation is given by 

(a) n (b) n x n 



2 F 1 (a,b; c;x) = ^ 



n=0 



nl 



(3.23) 



where (*) n = T(* + n)/T(*) is the Pochhammer symbol. In order to make a contact with the 
normalization used by Klemm et al. ,@ we rescale as 



CJi = 2y/uAUi, 0J 2 = U>2, 

U 

~ vA ~ 

= 2V«M2i, fi 2 = — ^2- 
u 



(3.24) 



Then the periods a,- and an. can be given by 



di = - (uj\ + uj 2 ) , a 2 = - (uJi — uj 2 ) 



a Dl = -- (fix + 3fi 2 ) 



(5iUi — 8 2 u 2 ) 
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a D2 = —- {Q,x - 3fi 2 ) (<Wi + 8 2 uj 2 ) 



71 



(3.25) 



where 81 = i{5 — 3 In 3— 41n2)/4 and 8 2 = 3i(l+31n3)/4 are constants determined from asymptotic 
expansion of periods.! The identification of periods by our solutions can be easily established by 



expanding (|3.25|) at u = oo. On the other hand, for large v, i.e., on the patch P 2 , ipi are expanded 
at v = oo and then consider similar linear combinations. 

In this way, we can check that the series solutions for v = v\ and v 2 comprise in fact the four 
periods. 



D. The fifth solution 



We have seen that the four solutions with indicial indices v\ and v% of Burchnall's fifth order 
equation in fact yield the four periods of the SU(3) gauge theory. However, there exists an extra 
solution in ( |3.9| ). Therefore, the appearance can be regarded to be characteristic in the ordinary 
differential form of the partial differential system. 

In this subsection, we show that it is possible to derive a fourth order equation satisfied by this 
fifth solution with aid of the SU(3) Picard-Fuchs ODE and we derive the closed formula for this 
fifth solution as an application of the SU(3) Picard-Fuchs ODE. 

Firstly, let us rewrite (|3.9f) in the form 



d 5 



d 4 



dz* + + ' ' ' + CBfi 



0. 



(3.26) 



where 2 = 1/2 and cb,% are some functions in x, y and z. As we have already seen in the previous 
section, since the four of the five independent solutions to (|3.9| ) are regarded as the four periods 
of the SU(3) gauge theory, we can write the Wronskian for ( |3.26| ) as 

CL\ CL2 &Di ^Di h 
h' 



W, 



B 



a 1 a 2 a Dl a D2 



(3.27) 



aS 4) a?) ag 

where ' = d/dz and the fifth solution is denoted by h. Again from basic differential calculation, 
Wb is generated from cb,4 and is found to be (c.f ( f2.15p ) 

1 



W. 



B 



(3.28) 



zas/^[x 2 + (y-z) 2 -2x(y + z)] 3 
up to the normalization of integration constant, which is irrelevant to the following discussion. 

Next, recall that the fourth order derivatives of periods can be reduced to a linear combination 
of lower order derivatives by using the SU(3) Picard-Fuchs ODE. Therefore, from (|3.27|) with the 
SU(3) Picard-Fuchs ODE, we see that 



W B = W SU (3)(h^ + c 3 h"' + c 2 h" + Cl h' + coh), 



(3.29) 



where q are given by (|2.10|) . From ( 3.28 ) and ( [2. 17f ), it is immediate to obtain the differential 
equation for h 



+ c 3 h + c 2 h + c\h + coh = R, 



(3.30) 



where 



R 



1 



(3.31) 



z 37 / 12 [5z - 9{x - y)][x 2 + (y- z) 2 - 2x(y + z)\ ' 

It is interesting to note that h satisfies the SU(3) Picard-Fuchs ODE with a source term. It is now 
easy to get a general solution to (|3.30|) (see also appendix B) 
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where pi and 6j are integration constants, the integration symbol is the integration constant free 
integral and Wi are determinants defined by 



(3.33) 
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In (|3.32|) , we have written as Wsu(3) for simplicity, but actually it should be substituted by ( j2.i7| ). 

To summarize, we have succeeded to find a closed representation of the fifth solution by using the 
SU(3) Picard-Fuchs ODE. It is interesting to compare ( |3.32| ) with the fifth solution of other three- 
term differential equation discussed by Burchnall.il In the case of a product of two Bessel functions, 
for instance, the fifth solution is expressed by an integral of a product of two "Wronskians" each of 
which is written by other two independent solutions.il However, our fifth solution does not admit 
such factorization, so ( |3.32| ) seems to imply the fact that the Appell function F 4 with parameters 
a = (3 = —1/6,7 — 2/3 and 7' = 1/2 (or equivalently, the SU(3) periods) can not be factored into 
a product of two (non-trivial) functions. 

Remark: In the semi- classical regime, series representation of this fifth solution corresponds 
to v 3 = 19/12 in ( \3.1$ ), and this can be also seen from the terms not including pi and €i in ^3.3^ ). 



IV. SUMMARY 



In this paper, we have discussed the Picard-Fuchs equations appearing in N = 2 supersymmetric 
Yang-Mills theories in view of ordinary differential equation and realized Picard-Fuchs equations 
as a system of ODEs. This construction has given a new relation among periods, moduli and QCD 
parameter by using the Wronskian and this is a systematic way to get such non-linear relation 
among periods in higher rank gauge group cases. 

In the case of SU(3), we have also found that Burchnall's ordinary differential equation for 
Appell's F 4 is a candidate of Picard-Fuchs ODE by identifying the SU(3) QCD mass scale param- 
eter with the scaling variable used in Burchnall's observation and confirmed that the four of five 
solutions of Burchnall's equation in fact coincide with SU(3) periods. As for the fifth solution, 
it has been shown that it has a simple and closed form by using the SU(3) Picard-Fuchs ODE. 
Of course, even if we consider arbitrary Riemann surface of genus two and try to get a similar 
result for the fifth solution to Burchnall's equation, the derivation will be failed because we can 
not always have an equation in fourth order like SU(3) Picard-Fuchs ODE. Note that the SU(3) 
Seiberg-Witten curve is a specific choice and the appearance of AgutS) plays the central role in the 
discussion. 

Generalization of Burchnall's construction of ordinary differential equation from a set of partial 
differential equations to Picard-Fuchs equations in other gauge theories is straightforward, but 
we do not know whether Burchnall type equations exist or not for Picard-Fuchs equations in 
those gauge theories. Studying these cases will open further aspect of Picard-Fuchs equations 
and hypergeometric nature of the equations constructed from the homogeneous hypergeometric 
equations in multiple variables. 

Finally, as another direction, since our equation is ODE in contrast with the usual Picard-Fuchs 
systems, it may be possible to consider a relation to classical P^-algebras in view of Picard-Fuchs 
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APPENDIX A. DERIVATION OF THE BURCHNALL'S EQUATION 

In this appendix, we briefly review the derivation of the Burchnall's equation. The reader is 
also recommended to refer to the original paper.il 
Firstly, notice that from (|3.6|) it is easy to obtain 



+ 7 + 7' -2)9J y F = {xz9 y + yz9 x )XF, 



where X = (9~ + a)(9~ + (3). We can also obtain 



UF = [9~{9~ + 1 -l)-{x + y)zX]F 

= (2^ + 7-7')^, 
U'F= [6~(9~ + 7 ' -l)-( x + y)zX 

= (29 y + 7' -7) 9 X F 
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and 
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9 y {9~ + 7' - 1) (9~ + 7 - 1)F = [(9 7 + 7 - l)9 x 9 y + yz (9~ + 7 )X] F. 



Addition of the two equations in ( |A"3| ) provides 
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where X +1 = (0~ + a + 1)(6 I ~+ (3+1), and ( |AT| ) and (|A2|) have been used. Moreover, from (|A5|) 
we have 
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2(0~ + 7 + 7 ' -2)9 x 9 y F 



(A6) 



Therefore, the expected equation is given by rearrangement of 

(2#~+ 7 + 7 ' -2)WF = 2(#~ + 7 + 7 ' -2)VF, (A7) 



i.e. 



Y -xzY 1 X-yzY 2 X + {x-yYz\2e~+ 1 + 1 +2)XX +1 F = 0, (A 



where 



y = e~(e~+ 1 -i)(e~+ 1 - i)(0~+ 7 + 7 - 2)(2#~+ 7 + 7 ' -4), 

Yi = {0~ + 1 '){6~ + 1 + 1 ' -l){29~ + 1 + 1 ' -2) + 6~{6~ + 1 -l){29~ + 1 + 1 ')^ 

r 2 = (^ + 7 )(^ + 7 + 7 '_l)(2^ + 7 + 7 '-2) + ^(^ + 7 '-l)(2^ + 7 + 7 '). (A9) 



Of course (|A8|) is equivalent to ( p!9| 



APPENDIX B. GENERAL SOLUTION TO FOURTH ORDER ODE 

This appendix reviews a construction of a general solution to the fourth order linear ordinary 
differential equation 

+ P{x)y"' + Q(x)y" + R(x)y' + S(x)y = T(x), (Bl) 

where ' = d/dx and P,Q,R, S and T are some functions in x. 

Firstly, let T = and let y^ (i — 1, • • • , 4) be the fundamental solutions to 

y {4) + P(x)y"' + Q(x)y" + R(x)y' + S(x)y = 0. (B2) 

Then we assume that the general solution to ( [BID is represented in the form 

y = j2C i (x)y l (B3) 
i=i 

by using unknown coefficients Cj. Differentiating (|B3|), we can obtain 

V =jZC' l y l + Y J C l yi (B4) 

i=l i=l 

but we further assume that the first term in the right hand side vanishes. Namely, we have 



E C'iVi = (B5) 



and 
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Repeating differentiation and imposing the vanishing of terms including C i: we get 
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Then from (|B9|) and 0, we get 
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In this way, we can arrive at the matrix equation determining all Cj 
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Consequently, Cj are given by 



(B12) 



d = a - ( -rrh d ^i C 2 = c 2 + f -r^-dx, C 3 = c 3 - ( -rzhdx, C 4 = c 4 + f —^-dx, (B13) 
Jo det r Jo det Y Jo det r jo det Y 

where q are integration constants, the integration symbol is the integration constant free integral 
and 
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